Scaling law for out-of-equilibrium smoothening of rough vicinal steps. Application to 

oxide surfaces during annealing 
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Langevin formalism approach is used to derive analytically the time evolution of the correlation 
function of steps in vicinal surfaces during annealing back to thermal equilibrium. A new generic 
scaling form is obtained for smoothening far-from-equilibrium one dimension objects. This approach 
allows interpreting a series of AFM images of sapphire surfaces showing the thermal evolution of 
initially rough step edges. From this analyse, we deduce that the step smoothening for (1,1,0,2) 
vicinal surfaces is governed by surface diffusion via a detachment-attachment mechanism during 
annealing at 1273K (~ O.GT™ 

citing ) • 
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Understanding how the competition between elasticity 
and disorder sets the morphological dynamics of inter- 
faces is a central theme of current research in statistical 
physics. It appears in the context of surface growth [1, 2], 
thermal fluctuations of vicinal surfaces [-3], domain wall 
motion in magnetic [4] or ferroelectric [-5] materials, flux 
lines in superconductors [6], cracks in solids [7], among 
other realizations. One of the most striking results in this 
field is the observation of spatio-temporal scale-invariant 
features characterized by exponents remarkably similar 
in systems with no apparent connection between them 
(see [1] for a review). Guided by this universality, it was 
then suggested that the interface motion should be de- 
scribed by the simplest possible Langevin equations com- 
patible with the symmetries of the problem [8] . 

This Langevin formalism was successfully applied to 
describe the fiuctuations of crystal surfaces at thermal 
equilibrium. This approach was used to interpret the spa- 
tial and temporal correlation functions measured from 
STM images [-3, 9, 10, 11] and allowed to measure quan- 
tities as e.g. the elementary kink energy [9, 10] or atomic 
diffusion energies [3, 11]. This is central for understand- 
ing and predicting the stability of surfaces/interfaces, 
and therefore the stability of associated nanostructure 
together with their physical and reactivity properties. 

Beyond thermal equilibrium studies, one may observe 
the evolution towards stability of initially fiat (low T 
equilibrated surface below its roughening transition Tr 
[12]) or rough surfaces (after initial preparation, high T 
annealing, growth or other surface processes). Such sit- 
uations are the most common and in the following we 
show that the Langevin formalism approach may be ex- 
tended to give interpretation of these very far from equi- 
librium states, leading to the determination of physical 
parameters (surface tension, relaxation time constant). 
For illustration of the method we applied this Langevin 
formalism to the experimental study of a vicinal surface 
of oxides. Such surfaces are difficult to equilibrates as 
long annealing at high temperature are needed during 
which it is difficult to maintain the stoichiometry. In ad- 



dition, very few studies, except for scarcely AFM inves- 
tigations [1-3], concern the surfaces of insulating oxides. 
To follow the smoothening of the atomic steps, we derive 
analytically the time evolution of the spatial correlation 
function of the step position. This one is found to ex- 
hibit a generic form that generalizes the Family- Viseck 
scaling classically used to describe roughening of initially 
fiat interfaces. The predictions are then compared to ex- 
perimental results on vicinal (1, 1, 0, 2)-sapphire surfaces. 
The time evolution of the step height correlation function 
can be perfectly reproduced and values of both the step 
stiffness and the atom hopping rate are deduced. 

Langevin equation of the out-of-equilibrium interface 
dynamics - Specifically, we consider a vicinal surface such 
as the one depicted in Fig. 1. The vicinal angle being 
small, the step density is low allowing to neglect the 
step-step interaction and allowing the use of a 1-d model. 
The morphological evolution is characterized by the se- 
ries displacement h{m,y,t) of the step numbered m (for 
an ideally ordered system, x = mL, L being the nominal 
step-step distance) at time t and location y. Within the 
capillary wave model, the Hamiltonian H of an isolated 
step such as the ones represented on Fig. 1 can be written 
[14]: 

where rj refers to the stiffness parallel to the step edge. 
In the Langevin formalism, the step profile is assumed 
to relax locally to minimize its energy. The back force 
depends on the diffusion processes. Without matter con- 
servation, i.e. for evaporation-condensation of atoms at 
isolate steps, this force is simply proportional to the en- 
ergy gradient 5H/Sh{y), and the ID Langevin equation 
takes the form of an Edwards- Wilkinson equation [9, 10]: 

dh(y) r/T d^h , 

where F denotes the atom hopping rate and ^{y,t) is a 
noise term allowing for thermal fiuctuations. This noise is 
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FIG. 1: Sketch and notations introduced in the text: The 
frame {ex,ey} is chosen so that and Cy are perpendicular 
and parallel to the mean step edges respectively, if denotes 
the angle between terrace and the mean surface plane. 



uncorrelated in space and time and its correlation func- 
tion can be written: 



(3) 



Note that the very same constant F appears in Eqs. 2 
and 3 making the principle of energy partition at thermal 
equilibrium satisfied. 

To characterize the dynamics of a fluctuating in- 
terface h{y,t), one usually computes the time evolu- 
tion of the spatial correlation function G(Ay, t) = 
{{h{y + Ay,t) — h{y,t))'^y This function can be ana- 
lytically derived from Eq. 2. Let us flrst discretize 
the step proflle: hny{t) ~ h{y = nyay,t) {uy — 
{— iV/2, iV/2 — 1} with Uy the lattice constant) and 
call hg{t) = il/VN)J2„^K^it)e-'^''"i^ the y- discrete 
Fourier transform of h. fjq{t) is similarly deflned. Then 
Eq. 2 writes: 



dt 



Vqhq{t) +f}q{t) 



(4) 



with u — 7]T/kT. Its solution is given by: 

~ ——^ — Tf? t I ~ I % r, UQ U 

' (5) 

where hq{0) is the y- Fourier transform of the step proflle 
h{y, 0) at time t = 0. Back into real space, the step proflle 
can be computed, and the complete time evolution of the 
spatial correlation function G{Ay, t) is deduced: 

G{Ay,t) = 

2TayN ^ 1 A A 1 ^ 

" ^ ' 1 " > ' 1 - cos(27rAy — -) 
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an initially flat step under the action of the thermal noise 
^. This situation has been extensively studied in the con- 
text of interface roughening problems and G^°"^'^{Ay,t) 
was shown to take the speciflc form, so-called Family- 
Viseck form[Lj]: 



The form of G can now be discussed. It is written as the 
sum of two terms G = G™"9''(A?/, + G'""'°^''{Ay,t). 
The flrst term G™"9''(A?/, describes the roughening of 



Chough (^^y-f.-^ = t2/3/(Ay/tl/^) 

where f{u) = 



fo[u/uoY'^ ifu<Mo 
/o if U> Mo 



(7) 



Where the exponents C = 1/2, /? = 1/4 and z = C//^ = 2 
refer to the roughness, growth and dynamic exponents, 
respectively. The two parameters /o and uq can be related 
to the two parameters v and F of the Langevin equation 
(Eq. 2): /o ~ l.^Tv-^'^ and uq ~ l.lv'^/'^. The second 
term G""°''*''(Ay,i) of Eq. 6 describes the smoothen- 
ing of the initial proflle h{y,t = 0) because of the fine 
tension v in absence of thermal fluctuations. Without 
loss of generality, let us consider an initially rough un- 
correlated proflle with zero average and variance set to 
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then given by G 
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cos(27rA?/(g/ayiV))) which takes the new generic form 
(Note the minus sign of the exponents and the abscence 
of the ( exponent): 



G""°°*''(A?/,t) = t-2/3g(Ay/ti/^) 
where g{v) = 



go{v/vo)'^ ifv<wo 
go if v> vo 



(8) 



Where the two parameters go and wq can be related to 
the two parameters ly and F of the Langevin equation 
(Eq. 2): go ~ 0.2a^ayiy~^/^ and vo 2.91^^^"^. 

From the forms of G''°"9''(A?/, i) and G™°°*''(A?/, t), 
- given by Eq. 7 and Eq. 8 respectively -, one can de- 
duce the global behaviour of G{Ay,t). This is sketched 
in Fig. 2. Depending on the time t, two cases can be 
distinguished: 

• For small time, i.e. for t ^ aQay/T), the varia- 
tion of the spatial correlation function G{Ay, t) can 
be decomposed into three regimes. At small scales, 
thermal equilibrium is reached and G'"°"^''(Aj/, i) 
dominates G(Ay,i). At medium and large scales, 
G(Ay, t) results from the smoothening of the ini- 
tial roughness and G'""°°*''(Ay, t) dominates: 

Ay«A(t) G{Ay,t)^{r/,,)Ay 
X(t) < Ay < ^{t) G{Ay,t) « 0.025a^ay{i^t)-^/'^ Ay^ 
Ay » ^{t) G(Ay, t) « 0.2a^„ay{iytr^/^ (9) 



(6) 



with A(t) 



• For large time, i.e. for t S> (jQay/T, the influence 
of the initial conditions h(jj, t = 0) is not seen any- 
more. The spatial correlation function G(Ay, t) is 



FIG. 2: Schematic illustration of the spatial correlation func- 
tion describing the time evolution of a step profile starting 
from an initially uncorrelated rough profile h{y, t = 0) with 
zero average and (Tq variance, (a) Smoothing regime observed 
for small times (see Eq. 9). (b) Roughening regime observed 
at large times (see Eq. 11). 

then given by G"'°"f''(A?/, t) for all Ay with two 
power-law regimes: 

A2/«^(t) G(Ay,t)«(rHAy (10) 
^y » m G(Ay, t) « l.h{T /v){vtf'^ 

with £,{t) « {vtf/'^. 

Experiments - This can now be confronted to ob- 
servations performed on the smoothening dynamics of 
vicinal surfaces in oxide surfaces. The samples studied 
here consist in square slabs of sapphire (Le Rubis S.A., 
10mm side, 0.2mm thick) with a surface oriented close 
to the (1,1,0,2) plane (misorientation ~ 0.06°). They 
present important technological interest as substrate for 
nanostructures, magnetic thin layers or giant magne- 
toresistance devices [17]. These surfaces are first chemo- 
mechanically polished to an optical grade, and carefully 
cleaned in an ultrasonic bath. As a result, one gets a 
rough out-of-equilibrium vicinal surface. A series of an- 
nealing in air (essential to maintain the surface stoichiom- 
etry) at constant temperature for increasing cumulated 
durations are then performed in an oven. After each an- 
nealing phase, the sample is cooled down to room temper- 
ature (cooling rate of 30K/min) and the vicinal surface 
is imaged via a Molecular Imaging Pico+ AFM in con- 
tact mode with gold coated Si3N4 cantilevers (0.58 N/m 
stiffness) . Special attention was paid so that our setup al- 
lows for an accurate and reproducible positioning of the 
AFM tip onto the surface and allows to image the very 



FIG. 3: 5 X 5 /im^ contact mode AFM images of the topog- 
raphy of a (1102) vicinal surface of sapphire after a) 5400s, 
b) 12600s, c) 23400s et d) 37800s of annealing at 1273K. The 
axis coincides with the [1, 1, 0, 1] direction. Note the defect 
(black circles) in (a)-(c) that proves the accurate repositioning 
of the AFM tip between two successive annealing processes. 



same area after each annealing treatment (see e.g. [16] 
for details). 

Depending on the temperature T, two regimes can 
be observed. At low temperature 973K < T < 1173 K, 
the evolution of the surface morphology is governed by 
the coarsening of 2D islands through anisotropic Oswald 
ripening during the experimental time. This regime was 
studied in a previous paper [16] and will not be discussed 
thereafter. At high temperature 1173 K < T < 1473 K, 
all the islands shrink up and overhangs in steps disappear 
rapidly, in less than one hour for T = 1173 K. After this 
initial regime, the evolution of the surface morphology 
is governed by the smoothening of the vicinal steps that 
can be described by Eq. 2. 

Typical AFM images of the vicinal surface taken after 
increasing cumulated annealing durations are presented 
in Fig. 3. The meandering lines are the steps defining the 
terraces. Contrast enhancement, image analysis and edge 
detection allow measuring h{m,y,t) and then comput- 
ing G(Ay, t) = {{h{m, y + Ay, t) - h{m, y, t))^)„i,y where 
the average has been performed over all edges of a given 
image to improve the statistics. Note that the mean dis- 
tance between two successive steps remains constant and 
pretty large (~ 750 nm) . This suggests that the morpho- 
logical fluctuations of a given step is independent from 
the neighbouring ones, as assumed in the Langevin de- 
scription proposed in Eq. 2. 
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FIG. 4: spatial correlation function G{Ay, t) of the step edges 
of the (1102) alumina surface for various cumulative durations 
t of annealing at T = 1273/sr K. The axis are logarithmic. 
The symbols correspond to experiments while the continuous 
lines corresponds to the analytical expression given Eq. 6 with 
u = 7.2x 10"^ nm^s and T = 5.3 x 10"^ nm^s 



The experimental correlation functions are plotted in 
Fig. 4. They exhibit the scahng given by Eq. 8 (Inset of 
Fig. 4) expected in the far-from-equilibrium smoothen- 
ing regime. Initial conditions are set by the image taken 
at to = 5400 s, the step profiles of which give the initial 
Fourier amplitudes \hq\^[t[)). The subsequent experimen- 
tal correlation functions are then fitted successfully using 
the analytical expression (Eq. 6) (continuous lines in Fig. 
4) . This provides a severe test in favour of the present ex- 
tension of Langevin formalism to far out-of-equilibrium 
systems. 

Note that the linear regime expected at thermal equi- 
librium is not apparent after more than lOh anneal- 
ing. From the fitted values v = 7.2 x 10~^nm^.s~^ and 
r = 5.3 X lO^^nm'^.s^^ (Fig. 4), the initial roughness 
(Tq ~ 100 nm, and from the pixel size ay = 14.6 nm 
of the images, one estimates the time Teq to pass from 
the out-of-equilibrium smoothening regime to the equi- 
librium roughening regime: r^q = OQay/T ~ 32 days. 
This makes this last regime difficult to observe experi- 
mentally and emphasizes all the importance of the out- 
of-equilibirum extension of the Langevin formalism to de- 
termine the physical parameters in oxyde surfaces. 

In summary, we used here Langevin formalism to de- 
rive analytically the time evolution of the spatial cor- 
relation function of initially rough steps of vicinal sur- 
faces relaxing toward thermal equilibrium. Far from equi- 
librium, this one is found to exhibit a generic scaling 
(Eq. 8) independent of the details of the considered sys- 
tem provided that the dynamics remains linear. On the 
other hand, the value of the exponents {/3, z\ depends 
whether the dynamics and the noise are conservative or 



not (see [1] for review). In particular they are found 
to be {/3 = 1/8,0 = 4} for conservative dynamics and 
noise that is expected at the onset of active diffusion [3] 
and commonly observed in metalhc surfaces [11]. Such 
a regime is surprisingly not observed for the (1,1,0,2) 
sapphire surfaces investigated here. Two origins of these 
non-conservative processes can be invoked: (i) The chan- 
nelled structure of terraces which was shown to infiuence 
the Ostwald ripening at low T [16]; These channels par- 
allel to [1,1,0,1] and roughly perpendicular to the step 
edge makes the mobility along the step edge difficult. As 
a result, atoms emitted from the step are diffusing on the 
terraces then captured on the same step at a different site 
(detachment/attachment mechanism); (ii) The exchange 
of oxygen atoms, during anneahng, with the surrounding 
atmosphere, evidenced by a different morphology of ter- 
races after UHV annealing at 1273K [18]. This makes the 
extension of the present work to the various conservative 
cases and to higher dimensions of high interest. Work in 
this direction is under progress. 
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